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Abstract 

We calculate the transverse energy flow accompanying small x deep- 
inelastic events and compare with recent data obtained at HERA. In 
the central region between the current jet and the remnants of the 
proton we find that BFKL leading ln(l/x) dynamics gives a distinc- 
tively large transverse energy distribution, in approximate agreement 
with recent data. 



The structure function F 2 (x, Q 2 ) for deep-inelastic electron-proton scattering has recently 
been measured in the small x region accessible at HERA, x ~ 10 -3 . The measured values 
show a striking rise with decreasing x which is entirely consistent with perturbative QCD 
expectations based on the precocious onset of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
leading ln(l/a;) behaviour. However, the observed small x behaviour of F 2 can equally well be 
mimicked by conventional dynamics based on Altarelli-Parisi (GLAP) evolution ||, where the 
steep behaviour is either put into the starting distributions or, alternatively, generated by Q 2 
evolution from a very low scale Qq. 

To obtain a sensitive discriminator between BFKL and conventional GLAP dynamics we 
need to look into the properties of the final state. A relevant observable, which has been recently 
measured at HERA by the HI collaboration ||], is the transverse energy (Et) flow in deep- 
inelastic events. The deep-inelastic data hint at an excess of Et in the forward part of the central 
region when compared to Monte Carlo simulations which incorporate GLAP evolution. GLAP 
evolution corresponds to a summation of large In Q 2 terms, which is equivalent (in a physical 
gauge) to the summation of ladder diagrams with strongly ordered transverse momenta (h?) 
of the emitted partons along the ladder: that is Q 2 3> k 2 ^> ... ^> kf, where we have omitted 
the subscript T on k Ti . 

At small x it becomes necessary to resum the large ln(l/x) terms and this is accomplished 
via the BFKL equation. The gluon 'ladder' diagrams relevant to this equation do not have 
the strong-ordering of the transverse momenta that is present in the GLAP diagrams. As a 
result more transverse energy is expected in the central region (between the current jet and the 
proton remnants) than would occur from conventional GLAP dynamics. These expectations 
are confirmed by explicit calculations ||. They are also hinted at by Monte Carlo simulations 
(which incorporate BFKL effects) of the gluon radiation accompanying heavy quark production 
at sufficiently small x, see Fig. 8 of ref. || or Fig. 9(b) of ref. 0. 

In ref. [[J we were mainly concerned with formalism and in using analytic methods to gain 
insight into the general characteristic features of the BFKL description of the energy flow in the 
small x regime. For fixed a s we derived an analytic form of the Et flow in the central region 
of deep-inelastic events at small x. The E T distribution was found to be a broad Gaussian- 
shaped plateau as a function of rapidity, with a height that increases with decreasing x, and/or 
increasing Q 2 . We also performed numerical estimates of the Et flow, which included the effects 
of running a s . These, more physical, calculations qualitatively confirmed the characteristic 
features of the fixed-a^ treatment, but did not fully cover the central region. The BFKL gluon 
emissions are only one source of transverse energy. There will also be contributions arising from 
parton radiation from the current jet and the proton remnants, and some enhancement of the 
Et flow from the subsequent hadronization. Here, in order to attempt a realistic comparison 
with the recent data, we extend the BFKL-based calculations of Et to cover a larger interval 
of rapidity and we use a Monte Carlo to simulate the effects of radiation from the current jet 
and of hadronization. We choose the LEPTO Monte Carlo || , as it gives a good description 
of final state observables in deep-inelastic scattering (and e + e~ collisions) in regions insensitive 
to BFKL small x dynamics. The matrix element (ME) + parton showering (PS) structure of 
the programme has the advantage that the GLAP-based initial state radiation can be isolated 
and therefore, in principle, be substituted by BFKL gluon emissions. 

The energy flow accompanying deep-inelastic events, in a small interval about x, Q 2 , is given 
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where the transverse momenta are defined in Fig. 1(a). For simplicity we have omitted the 
longitudinal structure function and assumed that F 2 = 2xF±. It is straightforward to include 
the small correction arising from Fl = F 2 — 2xF\. The function / is the unintegrated gluon 
distribution of the proton in which the k 2 integration is unfolded. To be precise 

Xjg{x h n 2 ) = J -jjrf( x j, k 2 p , /i 2 ) (2) 

gives the conventional gluon density at a scale /i 2 . In the leading ln(l/xj) approximation / 
is found to be independent of the scale n 2 and for this reason we have omitted the scale 
variable from the arguments of / in (|l|). In this small Xj approximation the function / satisfies 
a BFKL equation which effectively sums the soft gluon emissions below the emitted (xj,kj) 
gluon in Fig. 1(a). The same remarks apply to T 2 such that for sufficiently small x/xj, the 
function T 2 becomes scale independent and satisfies a BFKL equation which effectively sums 
the soft gluon emissions above the emitted gluon. 

As in ref. [j|], f(xj, k 2 ) is determined for Xj < 10 -2 by step-by-step integration of the BFKL 
equation down in Xj starting from a gluon distribution at Xj = 10~ 2 obtained from the MRS 
set of partons of ref. |H|, and T 2 is calculated for x/xj < 1CT 1 , as in ref. jTl), by evolving down 



from the quark box (and crossed box) contribution, T 2 , evaluated at x/xj = 10 _1 . Even in 
the lowest x region accessible at HERA, x ~ 10~ 4 , the above QCD prediction is limited to 
the Xj interval 10~ 3 < Xj < 10~ 2 . In Fig. 2(a) the continuous curve for 10~ 3 < Xj < 10 -2 
shows the distribution for x = 10~ 4 and Q 2 = 10 GeV 2 - values which are representative 
of the lowest x regime accessible for deep-inelastic scattering at HERA. The infrared cut-off 
on the transverse momentum integrations is taken to be k^ = 1 GeV 2 throughout, a value for 
which the calculated [12| values of F 2 are consistent with the recent measurements at HERA 



||T|. The predictions for the Et flow are less sensitive to the choice of the cut-off than those for 
F 2 . The ultraviolet cut-off is chosen to be Q 2 / z as implied by energy-momentum conservation 
151, where z = x/xj. 



To extend the prediction of the Et flow into the region Xj > 10 we proceed as follows. 



First we continue to use formula (]!]), but with 
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with the gluon taken from ref. flQfl . This expression for / follows from (||]), since, in the leading 



ln(l/xj) limit, / is independent of /i 2 . However, as Xj increases we soon reach the stage when 
the scale dependence of f(xj, k 2 ,fi 2 ) can no longer be neglected and so (§) becomes invalid. The 
gluon distribution / should always be positive, whereas the logarithmic derivative of xg becomes 
negative with increasing Xj due to the increasing importance of the usual virtual corrections 
of the GLAP equation. For this reason (as well as the omission of a quark contribution) the 
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Et flow calculated from the above prescription will be a larger and larger underestimate as Xj 
increases above Xj ~ 10~ 2 . This effect can be clearly seen in Fig. 2. 

A better approach for these larger values of Xj is to use the usual strong ordering at the 
gluon emission vertex, fcj 3> k 2 , so that k 2 « fc|. Then, on making use of (0), eq.(0) simplifies 

to ' ' 

= ^/^^N^E/a^^W^i.^Q 2 ). ( 4 ) 

where the "effective" parton combination J2 a fa = <7 + §(<7 + ?) arises from the dominance of 
gluon exchange ||14|| . In this way we include the contributions when parton a of Fig. 1(b) is either 



a quark or an antiquark, as well as the gluonic component which was dominant in the smaller 
Xj region. For sufficiently large Xj (but away from the proton remnants) formula (|) will give 
a much more reliable prediction for the transverse energy flow, but as Xj decreases the strong- 
ordering assumption becomes less valid and the consequent neglect of regions of phase space 
means that the Et flow will again be underestimated. This is apparent from Fig. 2(a) which 
shows the predictions of (|J) as a continuous curve in the interval 0.01 < Xj < 1. In Fig. 2(b) 
we show the results for x = 5.7 x 10~ 4 and Q 2 = 15 GeV 2 . This choice represents the average 
values [15 1 of the variables for the Et distribution, accompanying the deep-inelastic events with 



x < 10 , which was observed by the HI collaboration |4|. Again there is a reasonably flat 
plateau in the central region, but about 0.4 GeV lower than that of Fig. 2(a) and covering 
a smaller rapidity interval on account of the larger value of x. In summary, in the region 
10~ 2 ~ Xj ~ 10" 1 the Et flow is underestimated by both (P and (|j), but at different ends of 
the interval. Formula (|I|) is valid for Xj ~ 10~ 2 , while formula (H) is reliable for Xj ~ 10 _1 . 
From the combination of the results shown in Fig. 2 we conclude that BFKL radiation (which 
accompanies deep-inelastic events) gives rise to an approximately flat Et distribution in the 
central region with a height, which increases slowly with decreasing x, of about 2 GeV per unit 
of rapidity, in the HERA small x regime. 

We also performed the above calculations using GLAP evolution along the ladders. Figure 2 
shows that, as expected, this evolution gives a much smaller transverse energy flow than BFKL 
evolution. The discontinuity in the GLAP results at Xj = 10~ 2 is due to the inclusion of the 
quark distributions in the calculation for the large Xj region. As expected the quarks have little 
importance at small Xj where the gluon dominates. 

It is useful to translate the energy flow dE?/ 'dlnxj into a distribution in terms of rapidity, 
y in the HERA frame. First we note that in the virtual photon-proton centre-of-mass (cm) 
frame the rapidity is given by 

2/(cm) = ^ ln (^)- (5) 

In regions where the distribution is reasonably flat, a good estimate of the y distribution is 
obtained if we insert the local average value of k 2 into (|5J). Secondly we translate from the 
virtual photon-proton cm frame to the HERA laboratory frame using the formula 

1 ( AxE 2 \ 

y-y(cm) ~ -In ( -qT) ( 6 ) 
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which is valid at small x, since then the frames are approximately collinear. E p is the energy 
of the incoming proton in the HERA frame. 

At this stage our numerical results, shown in Fig. 2(b), cannot be compared directly with 
the HI data [||]. The distributions of Fig. 2 correspond to Et arising from gluons radiated 
from the initial state. Admittedly, at small x, this is expected to be the dominant effect 
in the central region between the current jet and the proton remnants - the large energy 
flow that is predicted is characteristic of BFKL dynamics. However the calculation does not 
take into account the current jet and its associated radiation, nor does it include any effects 
of hadronization. We therefore need to estimate the importance of the various components 
contributing to the transverse energy flow. To this end we compare in Fig. 3 the HI measurement 
of the Et flow PJ with four different distributions obtained using the LEPTO Monte Carlo ||. 
The parameters of LEPTO have been tuned to the hadronization resulting from jet production 
in e + e~ collisions and, to some extent, tuned to EMC deep-inelastic data in the higher x regime 
[jj~6| . The four Et distributions resulting from LEPTO show the energy flow obtained (i) from 
only parton showers radiated from the current jet, (ii) with parton showers from the initial state 
incorporated, (iii) from only parton showers from the current jet but with hadronization effects 
included, and (iv) from parton showers from both the current jet and the initial state together 
with hadronization. We see that the Monte Carlo gives a good description of the current jet 
and its associated radiation, but that effects of initial state radiation and of hadronization are 
unable to give sufficient Et in the forward part of the central region at small x. This has 
prompted a further study of the LEPTO Monte Carlo in order to assess whether this deficiency 
at small x is genuine or if it can be tuned away [I6|. This requires looking in detail at the 



LEPTO modelling of the GLAP-based initial state parton showers, and of hadronization of 
proton remnants more complicated than those arising from a diquark. We note from ref. H 
that a Monte Carlo based on the colour dipole model appears to give a better description of the 
Et flow in the central region than does LEPTO. On the other hand the colour dipole model is 
less successful than LEPTO in describing energy-energy correlations 0] and the Q 2 dependence 



of the jet rates fl7|l . It remains to be seen how much is simply parameter tuning and how much 
is directly attributable to QCD dynamics in the various Monte Carlo simulations. 

In Fig. 4 we confront the BFKL predictions with the HERA measurements of the Et flow 
in the central region with x < 10~ 3 . The average values of the deep-inelastic variables for these 
data correspond to x — 5.7 x 10~ 4 and Q 2 = 15 GeV 2 ||15|| . The BFKL predictions are thus 



simply those of Fig. 2(b), but shown now in terms of y, the rapidity in the HERA frame. This 
translation is achieved via eqs ((5|) and (|6]). However, the comparison of the observed Et flow at 
small x with the BFKL-based estimates is clearly incomplete. As emphasized above, we have 
omitted the effects of hadronization and of radiation from the current jet. The magnitude of 
these effects can be estimated from the LEPTO Monte Carlo. The appropriate histogram of 
Fig. 3 is reproduced in Fig. 4. 

To obtain a first estimate of the total Et flow we could simply add the Et resulting from 
the BFKL emissions to the LEPTO distribution. In other words the hadronization effects, 
which in LEPTO arise from the colour string stretching from the current quark to the diquark 
proton remnants, are assumed to give an underlying rapidity plateau whose gross features are 
insensitive to the properties of the initial state radiation. The Monte Carlo results in Fig. 3 
give support for this assumption, since we see that the level of hadronization is approximately 
independent of whether or not we include parton showers from the initial state. 
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If such a straightforward addition were to be performed on the results shown in Fig. 4 then 
the total estimate of the flow would be in approximate agreement with the data. There are 
two caveats. First, there may be some danger of double counting of radiation from the current 
jet and, second, our simple additive treatment of hadronization may be too naive. However, 
the first problem does not effect the central region, and secondly hadronization effects appear 
to be much less than the BFKL signal. 

To conclude, we have shown that small x deep-inelastic scattering is accompanied by a 
large E? flow in the central region arising from soft gluon radiation. This is a hallmark of 
BFKL dynamics and arises from the relaxation of the strong-ordering of transverse momenta. 
The first experimental measurements of the Et flow in small x deep-inelastic events indicate 
that there is significantly more Et than is given by conventional QCD cascade models based 
on Altarelli-Parisi evolution. Instead we find that they are in much better agreement with 
estimates which incorporate BFKL evolution. The latter dynamics are characterised by an 
increase of the Et flow in the central region with decreasing x. Measurements of the energy 
flow in different intervals of x, in the small x regime, should therefore allow a definitive check 
of the applicability of BFKL dynamics. 
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Figure Captions 



Fig. 1: (a) Diagrammatic representation of formula ([I]) showing the gluon ladders which are 
resummed by the BFKL equations for / and JF 2 . (b) The representation of formula ([|), 
which is obtained by the simplification of (|]) when Xj is large; there is now strong-ordering 
at the (parton a)-gluon vertex. 

Fig. 2: Numerical calculation of the Et flow as a function of Xj for (a) x = 10~ 4 , Q 2 = 10 GeV 2 
and (b) x = 5.7 x 10~ 4 ,Q 2 = 15 GeV 2 . The latter choice of variables is relevant to the 
HERA data for x < 10~ 3 . The effects of the current jet (and its associated radiation) and 
of hadronization are not included. The continuous curves are based on BFKL dynamics: 
formula (fj) is used for Xj < 10 -2 and formula (|]) is used for Xj > 10 -2 . For comparison, 
the dashed curves show the Et flow calculated using GLAP evolution. The dotted curve 
shows the effect of including only gluons at large Xj in the GLAP evolution. 

Fig. 3: The data show the Et flow as a function of rapidity in the laboratory (HERA) frame which 
accompanies deep-inelastic events with x < 10~ 3 ||. The proton direction is to the right. 
The LEPTO Monte Carlo distributions correspond, in increasing order, to ME+PS (final), 
ME+PS(initial,final), ME+PS(final) + hadronization, and finally ME+PS (initial,final) + 
hadronization. 

Fig. 4: The data show the Et flow accompanying deep-inelastic events with x < 10~ 3 observed 
by the HI collaboration || in the central region. The continuous curves show the BFKL 
predictions of x = 5.7 x 10~ 4 and Q 2 = 15 GeV 2 , which correspond to the average values 
of the variables for the data sample. The histogram is the LEPTO Monte Carlo estimate 
from Fig. 3 of the effects of radiation from the current jet and of hadronization. 
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Fig. 1 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405400vl 
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